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A grand canonical system of hard-core bosons, subject to thermal fluctuations, 
is studied on a lattice. Starting from the slave-boson representation with fields 
for occupied and unoccupied sites, an effective field theory is derived in which a 
complex field corresponds with the order parameter of the condensate and a real field 
with the total density of bosons. Near the boundary between the normal and the 
superfluid phase we obtain the Ginzburg-Landau functional for the superfluid order 
parameter. A mean-field calculation shows that the critical temperature T c increases 
with increasing density up to a maximum and decreases with further increasing 
density. 
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I. INTRODUCTION 

Interacting bosonic quantum systems are of special interest because of the effect of Bose conden- 
sation. The classical example is the Bose-Einstein condensation in an ideal Bose gas @,§. The 
analogous phenomenon in a real (i.e. interacting) Bose gas is the transition from a normal to a 
superfluid state (e.g., in 4 He). Another related phenomenon is the condensation of a cold gas of 
bosonic atoms in a magnetic trap gj. There is a number of recent theoretical investigations of 
the effect of the interaction on the critical properties of the normal-superfluid transition based on 
Monte Carlo simulations ||[| and analytic calculations j^U- 

A fundamental model for the description of interacting bosons is the Ginzburg-Landau theory, 
usually called the Gross-Pitaevskii approach to Bose condensates ^Jl(J. However, there are dense 
systems of bosons (like the classical 4 He superfluid), for which the Ginzburg-Landau theory is valid 
only very close to the phase transition from the normal fluid to the superfluid, where the order 
parameter is small. Away from the phase transition the interaction of the order parameter field is 
more complicated than described by this theory. Moreover, the Ginzburg-Landau theory is only 
a description of the order parameter (which is related to the density of the superfluid) but does 
not take into account the interaction with the non-superfluid part of the system. Consequently, it 
does not provide reliable information, e.g., for the value of the critical temperature. The question 
of the effect of the interaction on the latter has been discussed in the recent literature in great 
detail and with controversal results ||. It seems that the critical temperature can be shifted by 
a variation of the density of bosons. In particular, at low density it was found that the shift 
AT c /To ~ (na 3 ) a , where To is the critical temperature of the ideal (non-interacting) Bose gas, n 
the total density of bosons and a the scattering radius of the hard-core interaction. Depending 
on the calculational method and approximations, the exponent a varies between 1/3 |Tl[] and 1/2 
p2] |. Recent Monte Carlo simulations ]||| support a — 1/3. This value was also obtained by a 
self-consistent calculation of the quasiparticle spectrum and in a 1 /N expansion || . At a high 
density the critical temperature reaches a maximum and decreases with even higher densities. The 
latter is a consequence of the depletion of the condensate due to interaction. 

In order to give a complete overview of the properties of an interacting Bose gas we need a 
model which takes fully into account all parts of the system of bosons. Close to the critical point, 
however, it should lead to the Ginzburg-Landau theory. Such a model was given by hard-core 
bosons, based on a slave-boson representation ]l3| ). Here we will briefly discuss this model and 
evaluate its critical temperature for different densities. 

The paper is organized as follows: In Sect. 2 the slave-boson respresentation of hard-core bosons 
on a lattice is introduced. Then in Sect. 3 two collective fields are defined. One represents the 
superfluid condensate the other the total density of bosons, as discussed in Sect. 4. In Sect. 5 
the total density, the condensate density and the critical temperature are calculated in mean-field 
approximation. Finally, concluding remarks and a discussion are given in Sect. 6. 
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II. THE MODEL: SLAVE-BOSON REPRESENTATION 



A continuous system of hard-core bosons with scattering length a is approximated by a lattice 
Bose gas with lattice constant a. Although this approximation is limited because it restricts 
configurations of bosons to be commensurate with the lattice structure, it is more suitable for 
the investigation of a dense system of bosons than the usually considered |$| 4 (Ginzburg-Landau) 
theory The representation of the model uses the slave-boson approach to hard-core bosons |lT 



(Originally the slave-boson approach was invented for the (fermionic) Hubbard model 14 Jl5f| . ) 
The latter relies on a picture in which a particle trades its position with an empty site on the 
lattice. Both, the particle as well as the empty site, are described by corresponding creation and 
annihilation operators. In a functional integral representation of a grand canonical ensemble of 
bosons with chemical potential \i at temperature T this can be formulated in terms of a complex 
field b x of bosons and a complex field e x of empty sites with the action 

S s .b. = ji^2K e xt x ,x'b x 'e x , - ^^mI^I 2 , 

x,x' X 

where the first term describes the exchange of bosons and empty sites in a hopping process at 
sites x and x' with rate t x<x >. A local constraint \e x \ 2 + \b x \ 2 = 1 takes care of the complementary 
character of the bosons and empty sites. In S s .b. we consider only the thermal fluctuations (i.e. a 
vanishing Matsubara frequency) because non-zero Matsubara frequencies are separated by a gap 
if the temperature T is non-zero. Here we assume that the temperature enters through the action 
Ss.b. but not through the constraint. Physical quantities can be calculated from the partition 
function 



Z = J e~ s -- b Y[5{\e x \ 2 + \b x \ 2 - l)dbxdb* x de x del 



For instance, the total density of bosons n is given as the response to a change of the chemical 
potential \i 

_ T dlogZ _ T 1 dZ 

where N is the number of lattice sites. To evaluate the density of the condensate p s we must 
add an external vector potential to the action S s .b. (i.e. a Peierls factor to the hopping matrix t) 
and measure the response to this potential We shall return to the corresponding expression 
subsequently. 



III. COLLECTIVE-FIELD REPRESENTATION 



Since the fields b x and e x are subject to the local constraint, we can not treat them in a con- 
ventional way as order parameter fields. It is necessary to eliminate the constraint in the partition 
function which can be achieved by integration over the fields. For this purpose we introduce a 
complex collective field <5> x and a real field (p x which break up the biquadratic term in the action 
S s .b.. This can be written as the (Hubbard-Stratonovich) transformation 

s,, 8 - r 5>(i - <)ZK + r £ . + £ ( J ) ■ ( *■ + _ *U ) (<). (2) 

xx ! x x \ / \ / \ / 

The new action S' gives the same partition function which can be seen by integrating over the 
collective fields. The field ip is necessary in order to invert the hopping term which is now the 
positive matrix 1 — t. Then we can perform the integration of the slave-boson fields. This becomes 
simple if the 2x2 Hermitean matrix in Eq. (||) is diagonalized by a unitary transformation. The 
latter leaves the constraint \b x \ 2 + \e x \ 2 = 1 invariant and gives the eigenvalues 

Ax,± =<Px + 1/2T - ,i/2T ± v/(^ + l/2T + M /2T) 2 + |$ x | 2 . 



Now the partition function reads 

Z = Z f e- Sb JJ(S(|& B | 2 + l e ^| 2 ~ l) e - T ^- x *-+^ 2 - x *<-^ 

J X 

xdb x db*de x de*d<p x d$ x d$>*, 

where the only non-local term is 

x,x f 

The e x and b x integration can be carried out (s. Appendix A) which yields 



Z = Z [ e- s »T\e- T ^^-^ - ,/,;, ,/-!', -/<!.; 

J *x,+ — A x ,_ 



Z o/ e " S6 Il e " 1/4T+M/2T / 



The constant Zo is the normalization factor of the Hubbard- Stratonovich transformation of Eq. 
(H). It is convenient to separate the field- independent factor ]J X 6 - 1 / 4T +m/2T tQ define the partition 
function 



* = Mil / + m/2T)» + 

i V(<P*+»m 2 + \**\ 2 

This partition function has the effective action for the collective field $ x 

<5 = Sb + Sq, 

where Sb is its non-local (i.e. hopping) part and 

So^iog^a^i 2 ) 



is its local (i.e. potential) part with 



dip. 



IV. INTERPRETATION OF THE FIELDS 

The introduction of the collective field has completely separated the hopping part Sb from the 
potential part Sq. The latter does not depend on the phase of the collective field. The hopping 
part alone describes free bosons (ideal Bose gas) with the usual complex field. This can be seen 
by writing the partition function of the ideal Bose gas as 

z ibg = n [ e e- n *< e *-"°>/ r ] = ni 1 - e-^-^/ T \-\ 

k rifc>0 k 

where k is a quantum number which characterizes the system of bosons. This can also be expressed 
in terms of an integral over a complex field as 

Zibg = I exp(-^[l- e -^^)/ T ]|$ fc | 2 )j]d<i> fc ^/7r. 

** i. i. 



(Notice that a rescaling of the field yields only a factor to Zjbg-) For the translational- invariant 
Bose gas k is the wavevector and we have the energy 

h 2 k 2 

£fe - — 

2m 

with the boson mass m. 

Using the Fourier components (1 — tfc) in the non-local term Sb we can compare the latter 
with the corresponding expression of the ideal Bose gas 

£[l_ e -(c»-M0)/T]| $fc |2_ (3) 
k 

By setting the hopping term (1 — t^) -1 equal to 1 — e~( 1 ' T '( €k ~ ,x °> we obtain 

t k = (1 - e (e*-W>)/T)-l. 

The chemical potential of the ideal Bose gas is restricted to po < whereas it can have any real 
value in the interacting Bose gas. In particular, we can choose po > such that tk > 0. Moreover, 
in the Bose gas we can now apply the continuum approximation 

V^E^-C'*^ / [-^(V^ + al^l 2 ]^, (4) 

where a — ^ x (l ~ \>- ^ * s important to notice that a = (1 — io) -1 is positive because 1 — t 
was defined as a positive matrix. This implies that Si, is a positive quadratic form. Sf, can be 
compared with the expression of the ideal Bose gas if we replace a by — /io- Then the expression 
in (Q) is an approximation of Eq. (|^) for small tk ~ Mch which applies, e.g., to high temperatures. 
This approximation is very common in the literature [ [l7|]Tg|l and can also be used in the case of 
hard-core bosons. 

Physical quantities can be expressed as expectation values of the new fields tp and <f>. For 
instance, we obtain from Eq. (|l|) for the total density of bosons the expression (see Appendix B) 

n = -+T(tp x ), 

where 



Thus ip x , which is conjugate to \e x \ 2 according to Eq. @, is related to the total density of bosons. 
Conversely, & x , which is conjugate to e x b* x according to Eq. (||), corresponds to the density of the 
condensate |l3|,[l6| and can be expressed as the expectation value 



2 , 



We notice that the potential part of the partition function So is symmetric with respect to /J, — > —/j,. 
This implies that the expectation value (<p x ) is an odd function of /i. Therefore, the density varies 
monotoneously with /i, as it should. The symmetry of Sq with respect to ji — > —fi implies that p s 
is an even function of p. This is a characteristic feature of our lattice hard-core bosons in which 
bosons and empty sites are dual to each other. 



A. Near the Critical Point 

The (/^-integration in Z\ can be performed numerically in order to obtain an effective potential 

-^log^d^l 2 ). 
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It is interesting to notice that the expression log(Zi(|$ :E | 2 ) is linear for large values of |$a;| (cf. 
Fig. 1). Therefore, the quadratic term in Sb suppresses the large fluctuations. This means that 
the \^> x \ 4 approximation of the Ginzburg-Landau theory provides a stronger suppression of large 
fluctuations than the complete theory. 
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FIG. 1. Behavior of the local part of the action S Ioc = aT\$\ 2 - log(Zi(|$| 2 )) for T = 0.1, a = 1.02 
and three different values fi — 0, 0.07, 0.1. There is a critical point at /i ~ 0. 



Near the critical point we can expand the free energy around \<fr x \ — because the order pa- 
rameter field $ is small. The result of this expansion is the Ginzburg-Landau functional for the 
collective field with 



where 
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(h\$ X \ 2 + h\$a 



(Ta -&i)|* x | a -&2|* a 
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(6) 
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Expression (^|) is also known as the Gross-Pitaevskii functional which is often used as a model 
of a dilute Bose system ||[L0"|, In our derivation we have /Jr and T-dependent coefficients. As 
one can see in Fig. 1 the term with |$| 2 can change its sign with a changing /x. On the other 
hand, away from <& = the behavior is controlled by the confining part |<I>| 4 . This is the typical 
Ginzburg-Landau picture of a second-order phase transition. 

The effect of the |<E>| 4 interaction on T c in the model of Eq. (p|) was recently investigated in detail 
in a self-consistent calculation [Q and in a l/N expansion [||. In this model it is assumed that 
the coefficient of the |$| 4 term depends on the the scattering length a and the total density by a 
factor a 3 n. This is valid at low density. The calculation shows that the critical temperature of the 
interacting system T c , normalized with the critical temperature of the ideal Bose gas, is shifted as 



1 - T /T c ~ c an 



1/3 



(7) 



TABLE I. The coefficients of the critical temperature shift in Eq. (M) from different works 
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Numerical investigations Q show that the critical temperature must descrease at higher densities. 
This depletion effect, which cannot be seen in the Gross-Pitaevskii approach, will be discussed in 
the slave-boson approach subsequently. 

V. MEAN-FIELD APPROXIMATION 

Since (p x is a field which appears only in local terms, it can be integrated out at each site 
independently for a given value of the condensate field <f> x . The treatment of <fr x is more difficult 
since it appears in non-local term Sb- It can be studied in terms of the classical field equation 

- T^V 2 + Ta- ' = 0, 

which is the extremum of the action Sb + So. A further simplification is the additional assumption 
that the condensate field varies only weakly in space: V 2 <I> w 0. This gives the mean-field or 
Thomas- Fermi approximation. Both densities, n and p Sl can be evaluated in mean-field approxi- 
mation. The mean-field free energy reads 

Fmf = -^ogZ = \ - + Ta|$| 2 - log^d*! 2 )), (8) 

where |$| must be at the minimum of Fmf- For a given chemical potential fi there is a critical 
value T c which separates two regimes: one regime for T < T c in which the minimum of the mean- 
field free energy is $| 2 > and another regime with |$| 2 = for T > T c . This can be seen by 
plotting Ta|$| 2 — log(Zi(|$|^)) (Fig. 1) which has, depending on fj,, either a minimum at |$| = 
or another one at |$| ^ 0. p3| The minimal value |$| varies continuously as one goes through T c a. 
The behavior of the densities as functions of /i at a fixed temperature is shown in Fig. 2. 
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FIG. 2. Mean-field result of the total density n (full curve) and the superfluid density p s (dashed curve) 
at Ta = 0.1. 



A. Near the Critical Point 

For a very dilute system the mean-field approximation is insufficient and the calculations of 
Refs. RM should be applied. However, at higher densities (more than n ~ 0.2) the mean-field 
approximation should be reliable. Then the critical temperature T c of the mean-field calculation 
is 

T c = bi/a = Z[(Q)/(aZ 1 (0)). 




The decreasing behavior of T c /Tq (To oc n 2 / 3 is the condensation temperature of the ideal Bosc 
gas) is shown in Fig. 3. Our mean-field result for T c /To agrees qualitatively with the Monte Carlo 
result of Ref. |§. However, we expect that fluctuations might reduce the critical temperature 
substantially close to n = 1. 
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FIG. 3. Mean-field result of the critical temperature T c of the hard-core bosons in arbitrary units, 
normalized with the critical temperature of the ideal bose gas To oc n 2//3 . This behavior is similar to the 
result found in a Monte Carlo simulation of Ref. [1] 

We can expand the mean-field free energy (^) in powers of |$| 2 up to I*!"! 4 . Then the minimum 
of F must satisfy the mean-field equation 



1*1 



T r -T 



T, - T 



TaZ^O) - Z[(0) _ 
TaZ'^Q) - Z'{(Q) ~ (T - T c )h - 2b 2 /a ~ ~2b 2 /a 



(T~T C 



(9) 



if the right-hand side is non-negative and |$| 2 = otherwise. Since b 2 < in a typical situation 
(cf. Fig. 1), there is a non-zero solution. The coefficient on the right-hand side of @ can be 
evaluated numerically for a given value of ji. 



VI. DISCUSSION 



The slave-boson representation is given by two fields b x , e. x which are subject to a constraint: one 
represents empty the other singly occupied sites. These two fields are replaced by two collective 
fields ip x and <& x which have a direct physical interpretation. The former couples to \e x \ 2 and 
the latter couples to the product e x b*. The introduction of the collective field & x enables us to 
integrate the slave-boson fields and to eliminate the constraint of the these fields. The effective 
field theory provides a "two- liquid" theory which is represented by the fields cp x and $> x . These 
fields correspond with the total density n and and the superfluid density p s , respectively. The 
latter is the order parameter of the condensation whereas the total density n does not indicate a 
critical behavior at the condensation point (cf. Fig. 2). Therefore, n can be fixed near the critical 
point and the theory can be expanded in terms of a small order parameter field <fr x . This yields the 
well-known Ginzburg-Landau (or Gross-Pitaevskii) theory with density-dependent parameters. In 
other words, there is a phase boundary in the \x-T phase diagram which separates the normal from 
the superfluid phase. In the vicinity of this phase boundary we can apply the Gross-Pitaevskii 
approach. 

The condensate field can be treated in mean-field approximation, assuming a homogeneous order 
parameter. This reveals the phenomenon of depletion of the condensate due to a strong interaction 
among the bosons which has been observed in superfluid 4 He. In the slave-boson representation it 
is a consequence of the duality of empty and singly occupied lattice sites which is reflected by the 
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constraint \e x \ 2 + \b x \ 2 = 1. Thus physical quantities are symmetric with respect to a half-filled 
system (i.e. n = 1/2 or /j, = 0). 

The lattice theory is not very accurate at densities n « 0.5 because it reduces the motion of 
bosons significantly in comparison with a continuous system. Moreover, the mean-field approxi- 
mation neglects vortices and flucuations of the order parameter. From this point of view we can 
only expect a qualitative agreement of our results with those, e.g., obtained from experiments. 
Nevertheless, the mean-field approximation should be reasonable if the density is not too low. The 
critical exponent of the order parameter should be renormalized due to fluctuations, using the 
renormalization group for the three-dimensional |$| 4 model. 

The main result is that the phase diagram of the slave-boson theory of the strongly interacting 
bose system has two normal phases (i.e. a dense and a dilute one) and a superfluid phase for 
intermediate densities. Near the transition points a Gross-Pitaevskii approach can be used to 
describe the physics of small order parameter field. The superfluid density is low in our hard-core 
system. This might be a consequence of the strong interaction which suppresses the superfluid 
density. 

In conclusion, we established an effective field theory which enables us to evaluate the properties 
of a strongly interacting system of bosons. It takes into account the order parameter of the conden- 
sate and the total density by interacting fields. It describes the phase transition between a normal 
phase and a condensed phase. The phase transition were studied in mean-field approximation. We 
evaluated the density-dependent critical temperature at densities n > 0.2 in which the mean-field 
approximation of the order parameter is reliable. 
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Appendix A 

The integration of the e- and the 6-field can be performed in Z for each point x independently 

J + M 2 - l)e-^^ 2 - x ^-^ 2 db x db* x de x de* x . (Al) 

The integrand does not depend on the phases of the field. Therefore, the phase integration con- 
tributes a factor An 2 . Moreover, we set s := \b x \ 2 and t := \e x \ 2 which yields for (A.l) 

OO />00 r- 1 

/ 5(s + t-l)e~ x ^^ x ^- s dsdt^TT 2 e -K,+t-K,-(i-t) dt 
Jo Jo 




Appendix B 

To write the total density of bosons we have to evaluate 

T dlogZ _ 1 T 1 dZ 
N dfi ~ 2 + NZ~djjl' 

Diffcrentation yields 

2 



Since 



dlM J ^ ' V(^ + ///2T)2 + |* 2 |2 



we can write for the previous expression 

^= / e •-- n - r » a smhW{ * + M/2T)2 + l<w **. - 

d ^ J V ¥ V(^ + M/2T)2 + |$ S | 2 

This implies 

z9m z^J** J V(^ + M/2T)2 + |$ £ p ^[11 Ul 3:1 'J 11 

= E^- 
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